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LONG-TERM  GOAL 

The  long-term  goal  of  this  research  project  is  to  improve  the  Navy’s  capability  to  provide  high- 
resolution  forecasts  of  phase-resolved  ocean  wave  and  current  fields  from  remotely-sensed  data  and 
our  understanding  of  subsurface  oceanic  motions  from  their  surface  manifestations. 

OBJECTIVES 

•  Develop  a  time-domain  model  for  simulating  highly  nonlinear  wave  interaction  with  flows  with 
an  arbitrary  distribution  of  vorticity 

•  Extend  the  model  to  include  the  effect  of  wave  breaking  and  wind-wave  energy  transfer  and 
validate  with  laboratory/field  data 

•  Extend  the  model  to  include  the  effect  of  density  stratification  and  nonlinear  surface  wave 
interaction  with  density  fronts 

APPROACH 

This  research  builds  on  VORTWAVE,  a  fully-nonlinear  wave  model  developed  by  Nwogu  (2009)  for 
describing  the  time-dependent  evolution  of  the  ocean  wave  and  velocity  fields  due  to: 

•  nonlinear  wave-wave  interactions 

•  nonlinear  wave  interaction  with  flows  with  an  arbitrary  distribution  of  vorticity 

•  wave  breaking 

•  wind-wave  growth 


VORTWAVE  solves  the  exact  kinematic  and  dynamic  free  surface  boundary  conditions  expressed  as  a 
set  of  evolution  equations  for  the  free  surface  elevation  and  tangential  velocities  at  the  free  surface.  A 
velocity-based  boundary  integral  method  is  used  to  close  the  system  of  equations  and  relate  the  normal 
velocity  of  the  free  surface  to  the  tangential  velocities.  The  computational  efficiency  of  the 
wave/current  model  is  significantly  improved  by  expanding  the  kernel  of  the  boundary  integral 
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operator  in  terms  of  a  wave  steepness  parameter  and  using  the  Fast  Fourier  Transform  (FFT)  technique 
to  evaluate  the  resulting  convolution  integrals. 

Given  that  VORTWAVE  solves  the  inviscid  equations  of  fluid  motion,  viscous  processes  such  as  wave 
breaking  and  wind-wave  energy  transfer  that  occur  at  scales  much  smaller  than  a  wavelength  need  to 
be  parameterized  in  tenns  of  the  large-scale  wave  motions.  The  effect  of  wave  breaking  is  simulated  in 
VORTWAVE  by  idealizing  the  whitecap  region  on  the  front  face  of  spilling  breaking  waves  as 
infinitely  thin  vortical  layers  that  exert  a  vortex  force  on  the  underlying  irrotational  wave  field.  The 
zero-shear-stress  free  surface  boundary  condition  is  then  used  to  relate  the  breaking-induced  vorticity 
to  the  kinematics  of  the  underlying  wave  field,  resulting  in  a  dissipative  pressure  that  is  proportional  to 
the  square  of  the  nonnal  velocity  of  the  free  surface. 

WORK  COMPLETED 

Our  initial  effort  focused  on  incorporating  a  Lagrangian  discrete  vortex  model  into  the  code  to  allow  us 
to  evaluate  the  surface  signature  of  subsurface  vortical  motions  in  remotely  sensed  images  of  the  sea 
surface.  Two-dimensional  numerical  simulations  were  initially  conducted  to  evaluate  the  ability  of 
VORTWAVE  to  reproduce  the  results  of  Telste  (1989)  who  investigated  the  nonlinear  interaction  of  a 
pair  of  counter-rotating  vortices  with  a  free  surface  in  initially  calm  water.  The  model  was  then  utilized 
to  investigate  the  interaction  of  a  vortex  pair  with  surface  gravity  waves  and  compared  to  the 
laboratory  results  of  Fish  (1991). 

We  next  worked  on  implementing  energy  dissipation  terms  into  the  code  to  mimic  the  effect  of  wave 
breaking.  The  FFT-accelerated  boundary  integral  wave  model  assumes  non-overturning  waves  and 
requires  an  explicit  specification  of  local  regions  in  space  and  time  with  active  breaking.  Data  from 
laboratory  experiments  of  Chiang  and  Hwung  (2007)  were  used  to  guide  the  development  of  a  criterion 
for  the  onset  of  breaking.  The  model  was  then  applied  to  the  post-breaking  evolution  of  spilling 
breakers  to  evaluate  the  ability  of  the  vortex-force  parameterization  to  describe  the  overall  energy 
dissipation  rate  and  energy  dissipation  process  as  a  function  of  frequency. 

RESULTS 

Nonlinear  Wave  Interaction  with  a  Pair  of  Counter-Rotating  Vortices 

We  present  sample  results  for  the  interaction  of  periodic  gravity  waves  of  wavelength  k  and  amplitude 
A  with  a  pair  or  counter-rotating  vortices  in  infinitely  deep  water.  A  computational  domain  of  width 
32k  was  set  up  with  a  Gaussian-shaped  wavemaker  used  to  generate  waves  inside  the  domain. 

Damping  layers  were  placed  at  both  ends  of  the  to  absorb  outgoing  waves.  Two  counter-rotating 
vortices  with  circulation  T  and  initial  separation  distrance  a  were  placed  at  an  initial  depth  h=5a  below 
the  still  water  level.  Figure  1  shows  a  2-D  map  of  the  stream  function  and  vortex  path  for  a  simulation 

with  Froude  number  Fr  =  Y  / yjga3  =  0.5  ,  wavelength  k=2 a  and  wave  steepness  kA  =  0.2,  where  k=2n/ 
k.  The  vortices  rise  at  constant  speed  wv  =  T  /  2n a  before  being  deflected  near  the  free  surface.  A 

close-up  view  of  the  velocity  field  and  free  surface  of  the  left  vortex  at  1 1  ( a 2  /  T)  =  47  is  shown  Figure 
2.  Snapshots  of  the  wave  profile  at  different  times  during  the  vortex  motion  are  shown  in  Figure  3.  The 
left  vortex  induces  a  strong  adverse  current  that  steepeens  the  waves  and  leads  to  the  “scaring”  of  the 
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free  surface.  The  wave  height  is  also  significantly  reduced  in  the  region  between  the  two  vortices. 
Simulations  at  higher  Froude  numbers  led  to  instabilities  due  to  wave  breaking. 
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Figure  1.  Vortex  trajectories  and  stream  function  at  t/(a2 ’/T)  =47  for  nonlinear  wave  interaction 
with  a  pair  of  counter-rotating  vortices  (Fr  =  0.5,  kA  =  0.2,  /Ja=2). 
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Figure  2.  Close-up  view  of  the  flow  field  at  t/(a2/r)  =47  for  vortex  pair  interacting  with  a  surface 

gravity  wave  (Fr  =  0.5,  kA  =  0.2,  )Ja=2). 
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Figure  3.  Modification  of  surface  wave  profile  by  a  submerged  vortex  pair 

(Fr  =  0.5,  kA  =  0.2,  )Ja=2). 


Onset  of  Wave  Breaking 

Chiang  and  Hwung  (2007)  conducted  an  extensive  series  of  experiments  on  the  evolution  and  breaking 
of  modulated  wave  trains  in  the  large  wave  flume.  The  incident  wave  train  consisted  of  a  carrier  wave 
with  amplitude  Ac  and  frequency  ooc ,  and  perturbations  at  sideband  frequencies  a>±  -  (oc  +  A  co  slightly 

detuned  from  the  carrier  wave  frequency.  Numerical  simulations  were  conducted  to  identify  a  wave 
breaking  criterion  to  implement  in  the  VORTWAVE  code  consistent  with  the  breaking  locations 
observed  in  the  laboratory  experiments.  Four  different  wave  breaking  criteria  were  evaluated  including 
the  local  wave  slope;  the  ratio  of  horizontal  fluid  velocity  at  crest  to  phase  speed;  the  ratio  of  horizontal 
velocity  gradient  at  crest  to  a  characteristic  frequency;  and  the  ratio  of  vertical  acceleration  at  crest  to 
gravitational  acceleration.  The  spatial  variation  of  the  maximum  values  of  the  different  breaking 
criteria  are  plotted  in  Figure  4  for  one  of  the  test  cases  with  initial  wave  steepness  of  s=C4c=0.15.  The 
numerical  model  did  not  become  unstable  at  the  experimental  breaking  location  as  would  occur  during 
overturning  in  mixed  Eulerian-Lagrangian  boundary  integral  models  but  further  downstream  when  the 
slope  of  the  free  surface  was  near-vertical.  The  local  wave  slope  |  rjx  \  at  breaking  is  fairly  close  to  the 
30°  angle  obtained  in  Stokes  limiting  wave  with  a  120°  corner  flow  at  the  crest.  The  kinematic  (u  / C) 

criterion,  however,  is  much  lower  than  the  often  suggested  value  of  1.0.  Compared  to  the  wave  slope 
and  relative  crest  velocity,  the  horizontal  velocity  gradient  and  vertical  acceleration  increase  more 
rapidly  just  prior  to  breaking. 
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Figure  4.  Spatial  variation  of  different  wave  breaking  criteria. 

Spectral  Signature  of  Wave  Breaking 

The  measured  water  surface  elevation  time  histories  in  the  pre-breaking,  active  breaking  and  post¬ 
breaking  regions  were  also  analyzed  to  gain  a  better  understanding  of  the  signature  of  wave  breaking  in 
spectral  space.  Since  different  power  laws  have  been  proposed  for  the  saturation  range  of  the  spectrum, 
the  spectral  densities  of  gauges  located  in  the  active  breaking  region  were  re-scaled  by  factors  of  of 
and  co 5  and  are  plotted  in  Figure  5  for  three  wave  steepnesses  (s  =0.13,  0.15,  0.17).  The  different 
spectra  appear  to  collapse  onto  a  universal  line  in  the  high-frequency  range.  While  the  co 4  -scaled  plot 
appears  to  show  an  co~4  for  co  >  3 coc ,  the  of  -scaled  plot  shows  a  subtle  transition  from  of 4  to  co'5  at 
co  «  5 coc .  The  results  provide  further  support  to  the  hypothesis  of  a  Kolmogorov-type  cascade  in 
breaking  waves  with  the  breaking  process  primarily  dissipating  energy  transferred  to  the  high 
wavenumber  modes  by  nonlinear  wave -wave  interactions. 


Figure  5.  Frequency-scaled  spectra  in  the  active  breaking  region  (a)  scaled  by  of  (b)  scaled  by  co5. 
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IMPACT/APPLICATIONS 


It  is  anticipated  that  this  research  project  will  lead  to  an  improved  operational  tool  for  the  real-time 
forecasts  of  ocean  wave/current  conditions  for  naval  operations,  techniques  to  identify  the  surface 
signature  of  submerged  vortical  flow  structures  and  better  parameterization  of  unresolved  processes  in 
phase-averaged  ocean  wave  models. 

RELATED  PROJECTS 

This  project  is  closely  related  to  the  “High-Resolution  Ocean  Modeling”  project  (N00014-12-C-0566) 
being  undertaken  by  David  Walker  at  SRI.  The  present  project  is  focused  on  developing  and 
incorporating  advanced  features  into  the  model  while  the  SRI  project  is  focused  more  on  assimilating 
remotely  sensed  and  other  in-situ  data  into  the  model  and  evaluating  the  model  performance  with  field 
data. 
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